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Abstract. A Diophantine m -tuple is a set Aof m positive integers such that afc + 1 is a perfect 
square for every pair a, b of distinct elements of A. We derive an asymptotic formula for the number 
of Diophantine quadruples whose elements are bounded by x. In doing so, we extend two existing 
tools in ways that might be of independent interest. The Erdos-Turan inequality bounds the discrep- 
ancy between the number of elements of a sequence that lie in a particular interval modulo 1 and the 
expected number; we establish a version of this inequality where the interval is allowed to vary. We 
also adapt an argument of Hooley on the equidistribution of solutions of polynomial congruences to 
handle reducible quadratic polynomials. 

1. Introduction 

A Diophantine m-tuple is a set Aoim positive integers such that a6 + 1 is a perfect square for 
every pair a, h of distinct elements of A. For example, the first Diophantine quadruple {1,3,8, 120} 
was found by Fermat. (The nomenclature refers to Diophantus, who found a set that has the anal- 
ogous property in the rational numbers, namely {1/16, 33/16, 17/4, 105/16}.) It was proved by 
Dujella [3J that there are no Diophantine sextuples and only finitely many Diophantine quintuples 
(while a folklore conjecture asserts that there are no Diophantine quintuples at all). However, there 
are infinitely many Diophantine m-tuples for 2 < m < 4, and thus it is an interesting problem to 
try to count how many there are beneath a given bound. 

Dujella showed [2J that the number of Diophantine pairs contained in [l,x] is asymptotic to 
X log X, while the number of Diophantine triples contained in [1, x] is asymptotic to -^x log x. 
He considered the same counting problem for Diophantine quadruples, obtaining for sufficiently 
large x the lower and upper bounds 0.1608x^/^ logx and 0.5354x^/^logx, respectively, for the 
number of Diophantine quadruples contained in [l,x]. The primary purpose of this paper is to 
establish the following asymptotic formula: 

Theorem 1.1. The number of Diophantine quadruples contained in [1, x] is given by the asymptotic 
formula 

Cx^'Hogx + 0{x^/^{\ogxf'^+^'\\og\ogxf"^), 
where C = 2^/V3r(|)3 ^ 0.338285. 

Note that the exponent of logx in the error term, | + ^ ~ 0.90237, is indeed slightly smaller 
than 1 . 

In the course of establishing Theorem II. 1[ we found ourselves needing to extend two existing 
tools from the literature to suit our needs; these extensions might be of interest in their own right. 
The first of these tools is the Erdos-Turan inequality, which gives a quantitative bound for the 
discrepancy between the number of elements of a sequence that lie in a particular interval modulo 1 
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and the expected number. We require a version of this inequality in which the target interval is 
allowed to vary. To set some notation, let u = {«„}, a = {an}, and P = {Pn} be sequences 
of real numbers; we are interested in counting how many elements m„ lie in the corresponding 
interval /?„] modulo 1. (So that these intervals modulo 1 are sensible, we make the restriction 
«n < /^n < «n + 1 for all n.) Define the counting function 

Zn = Zn{u; a, /5) = #{1 < n < at : m„ G [a„, (mod 1)} (1) 

and the discrepancy between the counting function and the expected number 

N 

Dn = Dn{u; a, (3) = Zn - - a„). (2) 



n=l 



For any sequence {sn}, let Vn{s) = X]n=i W^n+i — Sn\\ denote its total variation considered as 
a sequence modulo 1, where \\y\\ = min^g^ \y — z\ is the natural metric on M/Z. We prove the 
following "moving target" extension of the Erdos-Turan inequality in Section [2l 

Theorem 1.2. Let (^^^d {Pn} be sequences of real numbers with an < Pn ^ c^n + 1- 

With Djq = Di\/{u; a, P) as defined as in equation we have 

\^n\ < + E (1 + ^KVN{a) + Vj,m) + -J M^{h) (3) 

h=l 

for any positive integers N and H, where c = 16/77r and 



Mfjih) = max 

l<T<Af 



T 



^^eihun] 



n=l 



(4) 



The following corollary is an immediate consequence of Theorem 11.21 and is easy to apply in many 
concrete situations. 

Corollary 1.3. Let {«„}, {a„}, and {Pn} be sequences of real numbers with a„ < /5n < «„ + 1. 
Suppose that both a and P are monotone sequences. With = Dn{u; a,P) as defined as in 
equation (|2]), we have 

N " 
\Dn\ ^Jj + (1 + |aiv-ai| + \Pn- Pi\)J2^n(^) 

h=l 

for any positive integers N and H, where M^^h) is defined in equation dH). 

We compare our Theorem 1 1.21 to the standard Erdos-Turan inequality, and exhibit two examples 
that probe the sharpness of our inequality, at the end of Section [2l 

The second tool that we extend is a result of Hooley on the equidistribution of the roots of 
polynomial congruences. Specifically, given any polynomial /(t) G Z[t], we consider the sequence 
of real numbers |, where k runs over all positive integers up to some bound y and u runs, for each 
k, over the roots of the congruence /(z/) = (mod k). Corresponding to such a sequence, we 
define the exponential sum 

Rf{h,y) = Y, e^-^'^^IK (5) 

^<y fM=0 (mod k) 
Q<u<k 
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Hooley [l4l[5l established nontrivial upper bounds for Rf{h, y) when / is an irreducible polynomial 
of degree at least 2. By adapting the methods of [4], we prove the following analogous upper bound 
for reducible quadratic polynomials in Section [31 

Theorem 1.4. Let f be a reducible, nonsquare quadratic polynomial with integer coefficients, and 
let D be its discriminant. Let Rf{h, y) be defined as in equation (|5]). For all real numbers y > 3 
and for every integer h 0, 



Rf{h,y) « f 1 + yilogy)^'~\\oglogy) 

P\h ^ 



5/2 



The trivial bound for Rf{h,x), namely the number of summands, has order of magnitude y log y 
(see Lemma 1531) . and so Theorem 11.41 represents a nontrivial upper bound for the sum when y is 
sufficiently large. It follows immediately from Weyl's criterion [7, page 1] that the normalized 
roots u/k are equidistributed modulo 1 . We remark that for reducible quadratic polynomials /, the 
true order of magnitude of Rf{h, y) is probably y\ therefore the estimate in Theorem II .41 cannot be 
improved too much. We discuss these issues further in Section |3l 

In Section |4] we show how these two tools can be used to prove Theorem 11.11 Several times 
in Sections |3] and m we need to invoke standard results on sums of multiplicative functions; to 
preserve the flow of ideas, we defer the proofs of all such results to Section [S] 

2. Erdos-Turan with a moving target 

In this section we prove Theorem 1 1.2 1 We begin by recalling the standard approach to bounding 
the discrepancy using exponential sums and Selberg's "magic functions". We then derive a 
bound for the individual Fourier coefficients that arise when using Selberg's functions. At that 
point, we can use partial summation to finish the proof of Theorem 1 1.21 At the end of the section, 
we compare our "moving target" inequality to the standard Erdos-Turan inequality and exhibit two 
examples that probe the sharpness of our inequality. 

2.1. Bounding the discrepancy by exponential sums. We begin by defining, for any positive 
integer H, the trigonometric polynomial of degree H 

Bh(x) = V / f sin 2TThx + — ^ V f 1 e(hx), (6) 

where e{y) = e^^''^ as usual and 

fiy) = -il-y)cotny-^. (7) 
We have 5^(0) = 1/{2{H + 1)) and Bnih) =Oif\h\>H+l, while if 1 < h < H then 

BH(±h) = —^ -(l-J!^^tf(J!^]]. (8) 

^ 2{H+1)\ H+l^''\H+lJJ 

We then define, for any real numbers a and P satisfying a < P < a + 1, the two related trigono- 
metric polynomials 

I3;y) = (3-a + Bniy - P) + Bnia - y) 
Snioi, /3;y) = p - a - BniP - y) - Bniy - «), 
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which satisfy S^{0) = P - a ±1/{H + 1) and, for h ^ 0, 

S+{h) = BH{h)e{-h(3) + BH{-h)e{ha) 

= -Bj^(-h)e{h(3) - BH{h)e{-ha). (9) 

These trigonometric polynomials are useful one-sided approximants to the characteristic function 
/9; y) of the interval [a, [3] modulo 1, which equals 1 if there is some number z between a and 
[3 such that y = z (mod 1) and otherwise: 

Proposition 2.1. For any real numbers a and (3 satisfying a < (3 < a + 1, 

SH{a,l3;y) < x{a,f3;y) < S^{a,(3;y) 

for all real numbers y. 

Proof. This is the fundamental property of the "Selberg magic functions" S*^; see flVl chapter 1] 
for an exposition which uses the notation we have employed. □ 

Note that the definition ([B of can be written as = Yl,n<N x{(^n-, Pn', Un)- Following the 
approach to proving the standard Erdos-Turan inequality, we use Proposition l2.1l to write the upper 
bound 

Zn <^S^{an,Pn]Un) =^ (pn-Oin + J^^—+ ^ S^{h)e{hUn)j 
n<N n<N ^ l<\h\<H ^ 

(where we have singled out the constant term 5^(0) of 5^); the definition (|2l) of thus yields 

T + 

n<N \<\h\<H 

= jf^ + Yl Yl {BH{h)e{-hPn)+BH{-h)e{han))e{hun) 

n<N l<\h\<H 

by equation Interchanging the order of summation, we get 

Dn < Jf-TY + Y (bhOi) ^ e{hun)e{-hPn) + Bni-h) ^ e{hun)e{han] 

l<\h\<H ^ n<N n<N 

The same calculation using Sjj instead of yields the corresponding lower bound 

N 



(10) 



Dn> 



J- J2 (BH{-h)Y,e{hun)e{hPn) + BH{h)Y,e{hun)e{-han)\ (11) 

l<\h\<H ^ n<N n<N ^ 



H + 



If the sequences a and /3 were constant, as in the standard Erdos-Turan inequality, then we 
could now factor out the exponential sum J2n<N ^i^^n) and then bound the remaining sum over h 
at once. Instead, we bound the Fourier coefficients Bnih) individually and use partial summation 
to estimate the effect of the varying terms e{—hPn) and e{—han) on the exponential sums. 

4 



2.2. Bounding the Fourier coefficients. We begin by establishing an inequality between the 
cotangent function and a rational function that we will use to bound the function / from the previ- 
ous section. 

Lemma 2.2. For < y < 1, we have 

1 I m-y) 1 
TT cot Try H < ' 



l-y y 2 2-y 

Proof. We start with the classical equality [[II equation (4.3.91)] 

1^/1 1 

TT cot -KZ = h > 

z ^-^ \ z ^ n n 



1 1 

z z + 1 z — 1 — z"^ 



> 



1 1 

- 1 ^ - 

n=2 

1 ,^111 

2;^ > — = - + 

- 1 ^-^ — 1 z z + 



n=2 

11 1 „ A 1 11 1 3^ 



z z + 1 z — 1 ^ nr — \ z -z + l z — \ 2 

n=2 

where the inequality is valid for < 1. Substituting z = \ — y yields 

-TT cot Try = 7rcot(7r(l - y)) > h , 

l-y 2-y y 2 

which is equivalent to the statement of the lemma. □ 

Lemma 2.3. Let f be defined as in equation ([7]). For Q < y < 1, we have \ f{y)\ < f — l). 
where c 

Proof. It is equivalent to show that 

^\fiy)\ 



16 
77r- 



l-y 

for < y < 1. By Lemma [Z2l we have 

^\f{y)\ 8 / , 1 

TT cot Try 



< 



l-y 7y V i-yj 7y 

^ 1 ^ 3(1 - y) 1 8 _ 21y3 - 63y2 + 30y - 4 



.y 2 2-y) ly 14y(2 - y) 

The denominator of the right-hand side is obviously positive for < y < 1; it suffices to show 
that the numerator is negative in that range. The polynomial p{t) = 21t^ — 63t^ + 30t — 4 has 
negative discriminant —29,484, and so it has exactly one real root. Moreover, p(l) = —16 and 
limj^^oo piy) = CO, so that real root is greater than 1; in particular, p(y) < for all y < 1. □ 

We are now able to bound the Fourier coefficients Bnih). 

Lemma 2.4. Let be defined as in equation For h ^ 0, we have 

,^,| 1/2-c c 
\BHih)\<-{^^ + 

where c = 



4:\H + 1 \h\ 



Proof. From the formula ([8]) for the Fourier coefficients of Bh, we see that 

2 / I I \ 2n 1/2 



\BHm 



< 



< 



2{H+l) 
1 

2{H + l] 
1 



\h\ 



H+l 

\h\ 



f 



\h\ 



H+l 



H+l 



+ 



f 



\h\ 



H+l 



(12) 



^^c/l-\h\/{H + l) 



2{H+1)\ 2\ \h\/iH+l 



by Lemma [23l this inequality is equivalent to the statement of the lemma. 



□ 



2.3. Finishing the proof of Theorem 11.21 At this point, we need only to record the outcome of 
a partial summation argument to be fully prepared to prove Theorem 11.21 Recall that VAr(s) = 



'^n=i W^n+i — Sn\\ dcnotcs the total variation modulo 1 of the sequence s = {s„}. 
Lemma 2.5. For any sequences ^f^d. of real numbers and for any integer h, 



N 



y^^e{hun)e{-hsr. 



n=l 



< {l + 2TT\h\VN{s))MN{h), 



where M]\f{h) is defined in equation @i. 

Proof. Let Exih) = ^2^=1 ^{hun)- Using partial summation, we have 



N 



N 



n=l 



^e{hun)e{-hsn) = ^ {En{h) - En-i{h))e{-hsn) 

n=l 

N-1 

= EN{h)e{-hsN) - ^ En{h){e{-hsn+i) - e{-hsn)), 



n=l 



and so by the triangle inequality 

N-l N-l 

e{hun)e{-hsn) < \EN{h) \ + ^ \En{h) \ ■ |e(-/is„+i) - e{-hsn) 



n=l 



n=l 

N-l 



l + ^\e{-hs 

n=l 



^i) - ei-hsr. 



From elementary properties of the exponential function, 

\e{-hsn+i) - e{-hsn)\ = |l - e(/i(s„+i - s„))| = |l - e(/i||s„+i - 
furthermore, |1 — e{y)\ < 2'n\y \ by the mean value theorem. We conclude that 

N-l , N-l X 

^ e{hun)e{-hsn) < M^ih) I 1 + 2n\h\ ^ - s„|| j , 

n=l ^ n=l ^ 

which establishes the lemma. 



□ 



Proof of Theorem \T72\ Beginning with the upper bound (ITOI) onD^, 

N 



Dn< 



(l3H{h)Ye{hun)e{-hl3n) + I3H{-h)Ye{hun)e{han)\ 

l<\h\<H ^ n<N n<N ^ 



we use Lemma 1X41 to obtain 



N 



H + 1 



l<\h\<H 



4:\H+1 



n<N 



e{hun)e{-hj3n) + ^ e{hun)e{ha 



n<N 



where c = We now invoke Lemma [231 to see that 



N 



H+1 



l<\h\<H 



2 - c c 
+ ^ 



4 ViJ+ 1 



< 



N 



(1 + 27r|/i|VW(/3))M^(-/i) + (1 + 2rt\h\VN{a))MN{h) 
l + \ J2 (|-^ + ^)(2 + 27r|/.|(V^(a) + \/;v(/3)))M^(/.), 



l<\h\<H 

since M^^—h) = Miy{h). At this point the summands for h and — /i are equal, and so 

H 



N 



H + 



h=i ^ ^ 



This is precisely the upper bound claimed in the statement of Theorem |1.2i the lower bound is 
established by exactly the same proof, starting with the lower bound (fTTI) on D^. □ 

2.4. Probing the sharpness of Theorem 11.21 The standard Erdos-Turan inequality flVl equa- 
tion (23)] has the form 



\D 



N 



< 



' h=l ^ ^ n=l 



If we restrict a and /? to be constant sequences in Theorem [L2l the conclusion is that 



\D 



N\ 



< 



N 



H 



N 



h=l 
H 



H + 



h=l 



H + h 



c c 

+ 



H + 1 h J 1<T<N 



max 



Ye{hun) 



n=l 



In theory, this conclusion is weaker than the traditional Erdos-Turan inequality due to the presence 
of the maximum; in practice, however, the attainable bounds on the exponential sums that arise are 
increasing functions of T, and so nothing would be lost. (Our method also does not include the 
possibility of replacing the | with /5 — a, as the latter difference is not independent of n.) 

Moreover, in its full "moving target" generality, the term Mj^(h) in Theorem II .21 cannot be 
replaced by | J2n=i ^(hun) \ ■ To see this, we consider the sequence {«„} = {n/ (A^ + 1)} for some 
positive integer N, and we select the "obliging target" intervals bounded by the sequences 



{a4 = K-2-"} and = K + 2-"}. 
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(13) 



The total lengths of these intervals is bounded by 2, but the number of points in the sequence u 
that lie in the intervals [a,/3] is N; therefore > — 2 in this situation. However, the total 



variations of the sequences a and (3 are bounded by |, while the exponential sum J2n=i ^(^ 
'^n=i ^{hn/ {N + 1)) has the exact value of —1 for every integer 1 < h < N. Therefore, for any 
integer H between 1 and A^, 



Ur 



N ^ / 9 — r r\ " 

Y^{l + ^h{V^[a) + V^m)i--^ + l\ J2e{h 

h=l V + / n=l 



H + 1 



N 

[''Mr 

n=l 



H + l '\H+l h H 

h=i ^ ' 



+ - < - + ii". 



Since the right-hand side can be substantially smaller than N — 2 (when H is near v A^, for in- 
stance), it is impossible for the left-hand side to be an upper bound for D^. Theorem [L2] remains 
valid in this situation, though, since Mj\f(h) is approximately ^, whence the right-hand side of 
equation ^ has order of magnitude A^ log H. 

We can also consider the sequence = {n'^} for some real number < 7 < 1 and the same 
"obliging intervals" (fT3l) as before, so that the discrepancy is again at least N — 2. The terms on 
the right-hand side of equation ([3]) must therefore conspire to give a contribution whose order of 
magnitude is at least A^. The total variations of a and P are each A^''' + 0(1), while the exponential 
sum J2n=i ^{hrC) can be shown to be asymptotic to N^~^ e{hN^) / {27iih'y) , which has order of 
magnitude N^^^ / h. Therefore the right-hand side of equation ([3]) has order of magnitude 

h=l 

showing that the bound is both correct and reasonably tight in this case as well. 

Both of these examples provide extremely positive discrepancies Dn, but simply replacing the 
"obliging intervals" with their complements a„ + 1] yields extremely negative disrep- 

ancies instead. 

Different methods of bounding the Fourier coefficients Bnih) yield different constants, in equa- 
tion (|3]), in the numerators of the fractions whose denominators are H+1 and h. The constants 2 — c 
and c we have derived above compare favorably to the constants yielded by other methods. Still, 
we can immediately identify two ways one could improve these constants if such an improvement 
were desired. First, the constant | in Lemma [23] is approximately 0.363783; from a numerical 
study, the best possible constant in this part of the argument would be approximately 0.356113, 
a difference of about 2%. Even aside from this, both of the inequalities in equation (fT2l) can be 
noticeably improved for \h\ near H by elementary means. 

3. Equidistribution of roots of quadratic congruences 

Our goal here is to modify the argument in Hooley's paper to show that the roots of a reducible 
quadratic are equidistributed in the same sense as his paper. 

3.1. Number of roots of reducible quadratics. For any polynomial / with integer coefficients, 
define 

P/(m) = ^{1 < r < m: f{r) = (mod m)} 



to be the number of roots of / (mod m) in a block of m consecutive integers. By the Chinese 
remainder theorem, the function p/ is multiplicative for any /. We also use the common notation 
ordp(ri) to denote the multiplicity with which the prime p divides n, and we write \\ n when 

ordp(?2) = a. 

Recall that the content of a polynomial with integer coefficients is the greatest common divisor 
of its coefficients. A polynomial with integer coefficients is called primitive if its content equals 1. 
Our first lemma allows us to reduce the task of counting roots of polynomials (mod m) to the 
primitive case. 

Lemma 3.1. Given f{t) E Z[t], let W be the content of f and write g{t) = ^f{t) E Z[t]. Let p 
be a prime, and let 7 = ordpiW). Then for any positive integer a, 



P", if a < 7, 

p"'pg{p"-"'), ifa>j. 



Proof. The congruence /(r) = Wg{r) = (mod m) is equivalent (HI Theorem 2.3(1)] to g{r) = 

(mod (1^^)- When m = p" with a < 7, then {W, m) = and the congruence is equivalent to 
g(r) = (mod 1); this is satisfied by every integer, in particular by all p° of the integers between 

1 and p". When m = p°' with a > 7, then {W,m) = p^ and the congruence is equivalent to 
g{r) = (modp""'''); this is satisfied by pg{p^~'^) integers in every block of p"^'^ consecutive 
integers, in particular by p'^ Pg{p°'~^) of the integers between 1 and p". □ 

In our investigation of the values of Pg for primitive reducible quadratics g, we will need the 
following elementary but awkward result. 

Lemma 3.2. Fix integers a and 5. Suppose that y and z are integers satisfying either 

(i) y = z < 6 or 

(ii) mm{y, z} = 6. 

If a < 26, then the inequality y + z > a is equivalent to y > [^], while if a > 25, then the 
inequality y + z > a is equivalent to the conjunction 

mm{y, z} = 6 and max{y, z} > a — 6. 

Proof. We begin by assuming that a < 25. First suppose that y > [|] , which is equivalent to 
y > ^ since y is an integer. If (i) holds, then indeed y + z = 2y > a; while if (ii) holds, then 
y + z>y + 6>^ + ^ = a. On the other hand, suppose that y < [|] , which is equivalent 
to y < |. If (i) holds, then y + z = 2y < a; while (ii) is impossible due to the contradiction 
y < ^ < 5 = mm{y, z} < y. This establishes the first asserted equivalence. 

Now we assume that a > 26. The condition (i) forces both y + z < 26 < a and mm{y, z} < 
6, which makes both sides of the second asserted equivalence false. On the other hand, under 
condition (ii) the inequalities y + z > a and max{y, z} > a — 5 are equivalent, since y + z = 
rmn{y, z} + max{y, z}. This establishes the second asserted equivalence. □ 

Any primitive reducible quadratic can be written as g{t) = {at + b){ct + d), where the primitivity 
implies that (a, 6) = {c,d) = 1. Define A = \ad — bc\, and note that the discriminant of g 
equals A^. We note that g is the square of a linear polynomial if and only if A = 0. 
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Lemma 3.3. Let g(t) = {at + b){ct + d) where a,b,c,d E Z and (a, b) = (c, d) = 1, and define 
A = \ad — bc\. Assume that A 7^ 0. Let p be a prime, and let 5 = ordp( A). Then for any positive 
integer a, 



'2, 


ifp- 


' acA, 




0, 


ifP 


ac andp A, 




< 1, 


ifp 


ac andp\ A, 




pL"/2j ^ 


ifp- 


' ac and p A and a 


< 26, 


.2/, 


ifp- 


' ac and p A and a 


> 26. 



(14) 



In particular, Pg{p^) < 2p^. 

We remark that since (a, b) = (c, d) = 1, any prime that divides two of a, c, and A automatically 
divides the third. Therefore the second line of the formula (fT4l) is the case where p divides both of 
a and c, while the third line is the case where p divides exactly one of a and c. 

Proof. First suppose that p divides both a and c. Then the congruence g{r) = (mod p°) implies 
bd = (modp), which is a contradiction since {a,b) = {c,d) = 1. Therefore Pg{p°') = in 
this case. 

Next suppose that p | A. Since A is an integer linear combination of ar + 6 and cr + d for 
any integer r, it is impossible for p to divide both such integers simultaneously. Thus the number 
of roots of (y((r) = (mod p") is simply the sum of the numbers of roots of ar + 6 = (mod p") 
and cr + c? = (mod p^). The number of roots modulo of ar + 6 = (mod j9°) is 1 if p f a 
and if p | a (since p \ b m the second case), and similarly for the number of roots modulo 
of cr + (i = (modp"). Therefore Pg{p°') equals 1 if p divides exactly one of a and c and 2 if it 
divides neither. 

Having disposed of the first three cases, for the rest of the proof we may assume that p \ ac and 
p I A (so that 6 > 1). The congruence g(r) = (modp") is then equivalent to (r + ba~^){r + 
dc^^) = (modp"). We may translate r without affecting the number of roots modulo so it 
is equivalent to look at the congruence r(r + Ai) = (modp"), where Ai = dc^^ — ba'^ = 
±A/ac. To calculate pg{p°'), we thus need to count the number of integers 1 < r < p° such that 
ordp(r) + ordp(r + Ai) > a. 

Note that p^ \\ Ai as well (meaning that p^ exactly divides its numerator while p does not 
divide its denominator). The p-adic ultrametric inequality tells us that of the three numbers 
{ordp(r), ordp(r + Ai), 6}, the two smallest are equal (or all three are equal). This implies that for 
any integer r, exactly one of the following situations holds: 

• ordp(r) = ordp(r + Ai) < 5; or 

• min{ordp(r), ordp(r + Ai)} = 6. 

If a < 26, we apply Lemma [X2l to see that ordp(r) + ordp(r + Ai) > a exactly when [|] < 
ordp(r). The number of such r between 1 and is just the number of multiples of pTo/si jj^^j 
range, which is exactly p'^ jp^°'l'^~\ = pL"/2J xhis settles the fourth case. 

On the other hand, if a < 26, we apply Lemma [X2l to see that min{ordp(r), ordp(r + Ai)} = 6 
and max{ordp(r), ordp(r + Ai)} > a — 6 > 6. For such integers r, either r or r + Ai must be 
a multiple of p"~^ (and not both, since a — 6 > 6). The number of such r between 1 and p° is 
therefore exactly 2p"/p"~^ = 2p^. This settles the fifth and final case. □ 

Lemma 3.4. Let f be a reducible, nonsquare quadratic polynomial with integer coefficients, and 
let D be its discriminant. Then p{m) < \fD ■ 2'^^™). 
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Proof. Note that D 7^ since g is not the square of a linear polynomial by assumption. As noted 
earlier in this section, we can write f{t) = W{at + b)(ct + d) for some integers W, a, b, c, d with 
(a, b) = (c, d) = 1. Let p be any prime, and let [3 = ordp(V^), which is an integer since D is a 
perfect square. Choose 7 and 5 such that p"' \\ W and p^ \\ {ad — be), so that 7 + 5 = /9. If a < 7, 
then Lemma [3n tells us that Pf{p°') = p°' < p"' < p^ . If a > 7, then Lemmas |3 . 1 1 and [33] tell us 

that = p"'pf{p°'~^) < p'^ ■ 2p^ = 2p^. In either case we see that p/(j9") < 2p°'''^p^'^\ Then, 

since p/ is multiplicative, 

p^(m) = W pfip"") < Y[ {2p^'M^)^j = Y[p"''^''^^'> < Jj2 = ■ 2^("^), 

p^ljm p\m. p\in p\m p\m 

as claimed. □ 

3.2. Proof of Theorem 11.41 and related remarks. We now describe our adaptation of Hooley's 
argument [4] to the case of reducible quadratic polynomials, which culminates in a proof of The- 
orem 11.41 After the proof we make some additional comments on possible improvements to the 
theorem. 

The only changes that need to be made to Hooley's argument Q involve the differences in the 
function pf resulting from the fact that / is now reducible. For one thing, the average size of p/ is 
now logarithmic rather than constant; this does not ruin the argument, but it needs to be taken into 
account. On the other hand, we know pf much more exactly for reducible quadratic polynomials 
(as Lemma [33] show s ) than we do for general irreducible polynomials; consequently, we are able 
to be more explicit in some stages. We also wish to make all dependencies on h and the polynomial 
/ explicit in our upper bounds, for the benefit of anyone wishing to employ Theorem 11.4! in the 
future. Even so, the number of changes is small relative to the several-page length of Hooley's 
argument. We have therefore, with apologies to the reader, decided not to include a self-contained 
proof but rather to indicate the necessary alterations to Hooley's proof. We will use some of the 
notation therein without definition when defining the notation is superfluous to the current account. 

Outline of proof of Theorem \1.4\ We begin by examining Hooley's auxiliary lemmata from [4J. 
His Lemma 1 -Lemma 3 are valid for any polynomial, reducible or irreducible, as the proofs es- 
sentially depend only on the Chinese remainder theorem. We will not use his Lemma 4-Lemma 
6: these deal with various estimates for the function pf, whereas we will simply use the informa- 
tion worked out earlier in this section. Finally, his Lemma 7-Lemma 8 make no reference to the 
polynomial and thus remain valid in our setting. (We note that since we are dealing always with 
quadratic polynomials, Hooley's parameters n and will equal 2 and 4, respectively, for us.) For 
the remainder of this proof, we use x instead of y so as to conform with the notation in [4j; this is 
not to be confused with the x that appears in the rest of this paper. 

Hooley's estimation proper of Rf{h, x) begins with a certain decomposition ^ equation (1)], 
namely Rf{h, x) = Si + S2, after which he observes that 

k<x 

We employ the upper bound in Lemma [J4[ to deduce that 

k<x 
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which differs from Hooley's estimate only in the presence of the v D term. Thus no further changes 
are needed in the estimation of S2: we obtain [4, equation (5)] 



x 

S2 < /d 



log x' 



which will turn out to be smaller than our estimate of Si. 

Subsequently, Hooley derives the upper bound (4^, equations (6) and (7)] 

In the estimation of S5, the only modification we need to make is to include a factor of D in the 
estimate 

(the factor of D resulting from applying Lemma [X4l as in the estimation of S2). Hooley's majoriza- 
tion of Se is exactly suitable for our purposes, except that we wish to keep explicit the dependence 
on h, so that we use his equation (9) rather than his equation (10). Our resulting version of JH 
equation (11)] is 

As Hooley does, we extend the range of summation on the right-hand side to all £ < x (the notation 
hides the fact that the sum currently runs over those integers less than x for which a certain divisor 
ki is at most x^/^). The resulting sum is treated in Lemma [53] below, whence we obtain 



Si < VD (logx)^ + — 

log X u V 



p\h 

This establishes the theorem. □ 

Theorem 1 1 .41 gives the estimate Rf{h, y) <^f^h Z/(log?/)^~^(loglog|/)^/^ for any nonzero inte- 
ger h. On the other hand, the number of terms in the exponential sum is J2m<y Pfi''^) ^fV l^Sl/ 
by equation (|29l) below. In other words, the exponential sum exhibits nontrivial cancellation for all 
nonzero h. By Weyl's criterion, this is precisely what is needed to show that the normalized roots 
of /(t) = (mod m) are equidistributed modulo 1. 

Note that for a linear polynomial f(t) = at + b with (a, b) = 1, the normalized roots are not 
equidistributed modulo 1, since they cluster around the points ^ with (s, a) = 1. A straightforward 
calculation and an invocation of Ramanujan's sum shows that 

k<y av+h=0 (mod k) ' 
0<i^<fc 

which is the same order of magnitude as the number of summands y, at least for some values of h. 
On the other hand, the roots of the reducible quadratic W{at + b){ct + d) modulo m certainly 
include the roots of at + 6 and ct + d, and so there will be a contribution to the exponential sum 
from these roots, whose order of magnitude is y. It seems reasonable to conjecture that the other 
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roots of the quadratic are distributed randomly. For example, with the polynomial /(t) = — 1, 
one would conjecture that 

k<y 1/2-1=0 (mod k) 
0<iy<k 

In any event, we should not expect any estimate of the form Rf(h,y) = o{y) for reducible quadratic 
polynomials /; therefore the bound in Theorem 1 1.41 is not too far from what one could prove. 

We remark that the dependence on D in the upper bound of Theorem 1 1 .41 could be improved if 
necessary. There are two places in the proof of Theorem 11.41 where we use Lemma [J4l to simply 
bound Pf{n) by \/13 ■ 2'^'^"^; but for many values of n the true size of Pf{n) is much closer to 2'^'^"^ 
A more precise application of Lemmas 13.11 and 13.31 would replace the \/15 in Theorem 11.41 by a 
smaller multiplicative function of D, one that was D^, for instance. Similarly, the dependence 
of the upper bound on h could be slightly reduced by using the exact expression derived on the last 
line of equation (|26l ) below. 

4. DiOPHANTINE QUADRUPLES 

With these two technical results in place, we can now proceed to the proof of Theorem 11.11 It 
turns out that the analysis hinges on studying a very specific family of Diophantine quadruples. A 
doubly regular Diophantine quadruple is one of the form 

{a,6,a + 6 + 2r,4r(a + r)(6 + r)}, (15) 

where a, 6, and r are positive integers satisfying a < b and ab + 1 = (so that {a, b} is a 
Diophantine pair); it is easy to verify that any such quadruple is in fact Diophantine. Let 

Q{x) = the number of doubly regular Diophantine quadruples contained in [1, x]. 

Dujella proved that almost all Diophantine quadruples are doubly regular; more precisely, he 
showed [2, Section 4] that the number of Diophantine quadruples contained in [l,x\ is Q{x) + 
(^^^0.292 |Qg2 Therefore it suffices to find an asymptotic formula for Q{x). 
Define the set 

R(m) = {y. I < u < m, = 1 (mod m)}. 

Notice that in any doubly regular Diophantine quadruple, we have the congruence r"^ = ab + 1 = 
1 (mod b). Moreover, the inequality a < b forces r < 6 as well, so that r must be an element 
of R{b). Conversely, any pair {r,b} with r G R{b) gives rise to a doubly regular Diophantine 
quadruple by taking a = (r^ — l)/6, except that r = 1 gives rise to a = which must be excluded. 
We must therefore count all the pairs {r, b}, with r e R{b) \ {1}, such that the largest element 
4r (a + r)(b + r) of the corresponding doubly regular Diophantine quadruple is at most x. In other 
words, if we define the function 

, s Jl, if 4r(a + r)(6 + r) < x, 
L(a, r, o; x) = < 

I 0, otherwise, 

then 

hm rei?(fe)\{l} ^ ^ 
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The obstacle we must overcome is this: whether or not 4r(a + r)(b + r) < x depends heavily 
on where in the interval [1, b] the congruential root r lies, when b is in the most significant range 
(around x^^'^ in size). By replacing the summand in equation (fT6l) by upper and lower bounds, 
Dujella was able [2, Theorem 3] to work out that the order of magnitude of Q{x), and hence 
of the number of Diophantine quadruples up to x, is x^^^ logx. We use our knowledge of the 
equidistribution of the roots r G R{b) to show that Q{x) is asymptotically the same as an analogous 
sum where the numbers r are chosen at random from [l,b]. 

We use the notation p(m) = ^R(m), which is a special case of the notation Pg{m) from the last 
section with g{t) = — 1. Applying Lemma [33] to this polynomial shows that 

{2, ifp^2orp" = 4, 
1, ifp" = 2, 
4, ifp = 2anda>3; 

in particular, 

p(m) = { 2^('")-i, if 21 II m, J> < 2^^'^)+^ (17) 




4.1. Truncating the infinite sum. To begin with, we need to truncate the infinite sum in equa- 
tion (fT6l ) in a manageable way. It turns out that we can accomplish this by sorting doubly regular 
Diophantine quadruples by a rather than b. 

Lemma 4.1. For any real numbers A,x > 2, the number of doubly regular Diophantine quadru- 
ples (fTSl ) satisfying a < A and 4r(a + r)(6 + r) < x is <^ (AxY^^ log A 

Proof. As before, the congruence = a6 + 1 = 1 (mod a) shows that r must be congruent to 
some u E R{a), so that we can write r = u + ak for some integer k. Furthermore, the inequality 
a < b forces a < r as well, so that k > 1. Conversely, any such r determines a doubly regular 
Diophantine quadruple by taking b = (r"^ — l)/a. Note that the inequality 4r(a + r)(b + r) < x 
implies that 

/(^y _j_ akV — 1 \ 
x>A{u + ak){a + u + ak)i ^u + akj > A{ak){ak){aP) = Aa^k^; 

in other words, it is necessary that 1 < A; < (x/Aa^Y^'^. We conclude that for every integer 
a and every u E R{a), there are at most (x/Aa^Y^^ corresponding doubly regular Diophantine 
quadruples contained in [l,x]. An upper bound for the number of such quadruples with a < A is 
therefore 

a<AueR{a) ^ ' a<A ^ ' a<A 

where the last step uses equation (l30l) . □ 

Let i\}{x) > 1 be any function that tends to infinity as a; — > oo, but more slowly than log x (we 
will choose a specific function i\){x) later in equation (|25l) ). The following proposition shows that 
we can truncate the sum in equation (fT6l) at 

B = B{x) = \_x^l^i){x)\. (18) 
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Proposition 4.2. For any real number x > 3, 

b<BreR{b) ^ ^ ^ 'TK J / 

where B is defined in equation (fT8l) . 

Proof. We begin by bounding the number of doubly regular Diophantine quadruples with b > B 

and 4r(a + r)(6 + r) < x. Since r'^ = ab + 1, these inequalities force x > 4r(a + r)(6 + r) > 
Ar^b > 4:{ab)b > AaB'^; in other words, every such Diophantine quadruple satisfies a < x/AB^ < 
x^f^ /il){x)'^ . By Lemma |4~T1 the number of such Diophantine quadruples is 



^■ip{xy J \ip{xy J ■ip{x)^/'^ 

Therefore equation (fT6l) becomes 



b<BreR(6)\{l} ^ ^ \ 7 / 



b<-BreR(fe) 

which is equivalent to the statement of the proposition. □ 
4.2. Invoking equidistribution. 

Lemma 4.3. Let a, b, r be positive integers with a < b and ab + 1 = r^. Then for any real number 
X >?>, the inequality 4r(a + r)(6 + r) < x is equivalent to ^ < \{b, x), where 

ifb < (x/16)i/^ 



Kb,x)={ 1x^2 I 1 b \ (19) 



263/2 4 



+ + ^6 > (a:/16)V3. 



Proof. The inequality 4r(a + r){b + r) < x is the same as 

If we set s = |, which is the normalized root used in the statement of our equidistribution result, 
then this becomes 

2/-I \2 - ^ s(l + s) X ^(b 



.^(l + .)^<^ + ^ = 5;;^l+0^-JJ. (20) 

Note that s < 1 and so the left-hand side is at most 4. If it happens that b < (x/16)^/^, then the 
first term on the right-hand side is at least 4, so that the inequality always holds; setting A(6, x) = 1 
is therefore valid in this region. Otherwise, for any positive real number y, we see by completing 
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the square that s^(l + sY < y if and only if s < (y^y + 1/4)^/^ — 1/2. Using this equivalence, the 
inequality (|20l) becomes 



X sil + s) 1 



463 4/ V Va;// 2 



1 \ 1/2 1 
X 1 \ ' 1 



7: + 



since b > 1. This establishes the lemma in the case 

wherefo > (x/16)^/3_ □ 

Proposition 4.4. For anj rea/ number a; > 3, 



where B is defined in equation (fTSl ). 



Q(x) = ^p(6)A(6,a:) +0 a;^/3^/'(x)(loga;)^/2QQgiog3,)5/4^ 



1 1/2 



Proof. Consider the concatenation of B finite sequences, the 6th of which consists of the roots of 
= 1 (mod b) normalized by dividing by b\ in other words, consider the sequence i?(l) U |-R(2) U 
|i?(3) U ■ ■ ■ U ^R{B). This sequence has S{B) elements, where we define 

b<y 

We will apply the moving-target Erdos-Turan inequality, Theorem |1.2l to this sequence; our target 
intervals will be [0, A(6, x)] for each element r/b of lR{b). With this setup, we have 

b<B b<Br€R{b) ^ ^ 

by Lemma |431 whence 

^5(B) = $^ J2 L('^,r,b;x) -Y,Pimib,x) -0). 

b<Br(^R{b) ^ ^ b<B 

In other words, 

g(x) = J] A(6, x) + O (d^cb) + (V^(x) + 



ip{x) 



by Proposition |4]2] (recalling that B is defined in equation (fTSi)). 

Our intervals [0, A(6, x)] have the property that the lower and upper endpoints each form mono- 
tone sequences; the variation in the lower endpoint is 0, while the variation in the upper endpoint 
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is 1 — A(i?, x) < 1. By Corollary II .31 we therefore have 

DsiB) « ^ + (1 + + 1) Ms^B){h) « ^ + E MsiB){h) 



h=l h=l 



for any positive integer H; we will choose H later in (1221) . By Lemma [531 the first term satisfies 

S{B) BlogB x^/^ij{x)\ogx 
H ^ H ^ H ■ 

By Theorem O applied with ?/ < 5 and with f{t) =t^ - I, so that D = 2, 

MsiB){h) « n + ^) ■ B{\ogB)^'~\log\ogBf' + max p(&), 

p\h V / 

the final term arising because we have to consider what happens if we chop off the exponential sum 
in the middle of one of the finite sequences lR{b); here even the crude bound p(b) < B suffices. 
Consequently, by Lemma [541 we have 

5^ MsiB){h) « HB{\ogB)^'-\\og\ogBf\ 

h<H 

We conclude that 

Q(a;) = ^p(6)A(6,a;) + o( '^'^'^^^^^"^'^ +gi?(logi?)^-^(loglogi?)^/^ + a:^/^(^(x) + 

b<B ^ 



V'(X)V2 

gp(6)A(6,a:) + o (^^^^^^^^|^ + i7a:V3^(x)(loga;)v^-^(log logx)^/^ + ^J^) • 



(21) 



To optimize this error term, we choose H so that the first two terms are the same size. This 
choice turns out to be 

H = [(log x)^-^/'(log log x)-^/^] , (22) 

and with this choice, we get 

Qix) = J2pib)X{b,x) +o(x'/'ij{xK\ogx)^'/\\oghgxf' + ^J^^^ 



b<B 



□ 



4.3. Calculating the weighted sum. We evaluate the sum 

5^p(6)A(6,x) 

b<B 

by partial summation, recalling that B = [x^/^^(x)\ where ^(x) tends to infinity but more slowly 
than logx; we also recall the notation S{y) = J2b<yPib)- After working hard to evaluate the 
eventual leading constant in closed form, we finally finish the proof of Theorem 11.11 at the end of 
this section. 
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Proposition 4.5. For any real number x > 3, 



b<B 



:x 



Proof. For notational convenience in this proof, we set Bi = Using the definition ([T£ 

of A (6, x), we have 



b<B b<Bi Bi<b<B 



Bi<b<B ^ ^ 

= E P(b)^^ + \ + P(B,)-^^SiB) + 0{l). (23) 

Bi<b<B 

The remaining sum can be written as a Riemann-Stieltjes integral, to which integration by parts 
can be applied: 



, ^1/2 1 pB / 1/2 1 \ 1/2 

Bi<b<B JBi \ / 



a;i/2 / a;i/2 i\ 1/2 



We note that 



V2B3/2 4y V2^(x)3/2 '4; 2 \^{xy/\ 
while the similar term with B replaced by Bi evaluates to exactly |; we also note that 

d_ fx^ 1\ _ _3x^ / 1\ 
rft V2t^/' 4^ ~ 8t5/2 V2t3/2 + 4 J 
Therefore equation (l24l) becomes 

E + i = (^O(^(x)-V^))S(B)-iS(i30+|-"' rUflf^+i)""*' 

Bi<b<B -^^1 ^ ^ 

which results in a lot of cancellation in equation (l23l) : 

cr+^ / ^1/2 1 \ -1/2 



6<B 

_ 3 



^BlogB 

il)[x)'^/'^ 
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by Lemma [531 Finally, we extend the integral to infinity, noting that 

-1/2 



so that 

J2 pmib, x) = \x'i^ j^^ ^ + 1 j dt+oi^ 



b<B 



V'(x)3/2 ^1/2 



Since B = [x^^^ip{x)\, both error terms are <^ x^/'^(logx)/^/'(x)^/^, which establishes the propo- 
sition. □ 

Lemma 4.6. For any real number x > 3, 



1 + ^1 "^^dt='^^x^/Hogx !\l + ^u)-'l^u-^l'du + 0{x'l') 



3x^ 

4 -/(./i6)i/3 ' f^^' J t^^' Jo 



Proof. First we use the asymptotic formula for S{t) in Lemma [531 



(z/16)i/3 



1 + ^t:^ ^ dt 



t3/2 J fo/2 



Since (1 + 2x^/'^ /t^'"^)-^/'^ < 1, the last integral is bounded above by 2{{x/lQf/Y^/'^. Therefore 
3x1/2/-- / 2x^/^\-^'^ S{t) ^ 9x1/2/- / 2xV2\-i/2iogt 3 



We now make the change of variables u = x^^"^ /At^^"^ in the remaining integral, which yields 

9x^/2 r 2x^/2 \-^/' 



2tx 



roo / 9 l/2\-V^ 

/ 1 + ^ r'/'\ogtdt 

f\ x-l/2 4m , xl/3 / xl/3 



2'/'^'^V(ioga;) r(l + 8w)-i/V2/3rfM + 2 [\l + 8uY/'^-^du 



This establishes the lemma, on noting that the last integral converges to some finite constant. □ 
Lemma 4.7. We have 

o2/3 /-l o4/3 



/ (1 + 8u)-i/2m-2/3 ^ — ^ ^ 0.338285 
Jo ^ ^ 3r(|)3 



Proof. This lemma turns out to be an exercise in mining known results about special functions. We 
begin with Euler's integral representation [1, equation (15.3.1)] for the hypergeometric function 
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2^1, valid for all complex numbers z other than real numbers greater than or equal to 1, as long as 

^c>m)> 0: 

b; c; z) = Jf^ f t'-\l - tr-'-\l - tzY^ dt 
T[h)T{c-b) Jq 

Choosing a = \,h = \,c= |, and z = —8, and recalling that T{y + I) = yT{y),v^c see that 

2^/^ r\. , 0^.-1/2^.-2/3.^. 2V3r(i)r(i) 3-22/3 



7r2 



^ (1 + 8.)-/2,-/3 = ''-/-^^.F. ih h I; -8) = '-^.F, (1, i; I; -S) . 



^3' 

Next, we apply the quadratic transformation [1 , equation (15.3.22)] 

2F1 {a,b;a + b + i; z) = 2F1 (2a, 2h;a + h+\-\- - z) 

to obtain 

3-22/3 , , 3-22/3 



7r2 >j o ' 



At this point, we can actually evaluate this special value by the formula [[11 equation (15.1.21)] 

2Fi(a, 6; a - 6 + 1; -1) = 2-V/2 ^ 

^ ' ' ' ^ r(l + f -6)r(| + f)' 

which gives 

^-IHfU 2.4. l^_ 3 r(|) _ 1 r(|) 

^2 2^1^1,3,3, 1; 2i/37r3/2r(|)r(l) 2V3vr3/2r(|)- 

We now invoke the duplication formula for the Gamma function [T^ equation (6.1.18)], namely 

T{2z) = (27r)-i/222--i/2r(2)r(^ + i). At z = | this yields r(|) = 2i/3vri/2r(|)/r(|), so that 

1 r(|) _ 1 r(|)2 



2i/3vr3/2r(|)r(l) 22/37r2 r(|) 

Finally, the reflection formula for the Gamma function (T, equation (6.1.17)] is T{z)T(l — z) 
Tc CSC Tcz; again choosing 2 = |, we obtain r(|) = tt csc |/r(|) = 27r/3"'^/2r(|), so that 

1 r(i)2 2^/3 



22/3vr2 r(|) 3r(|)3 

This chain of equalities establishes the lemma. □ 
Proof of Theorem [7Jl By Proposition 14.41 

Q(x) = Vp(6)A(6,x) + ofx^/3V^(x)(logx)^/2(loglogx)^/^ + ^^^^ 
By Proposition 14.51 

, 3X^2 .00 / 2x^/2 V'^' , ^/l/3,.^n ^./2/2n i ^5/4 X^/^logxA 

= ^ ^^^^^^ (l + ^j J^C^t+0(xV3^(x)(l0gx)v^/2(l0gl0gx)V^ + ^^j, 

At this point we choose 

tPix) = (loga;)(2-^)/3(loglogx)-^/^ (25) 
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to optimize the error term, yielding 
By Lemma l4!6l 

92/3 rl 

Q(a;) = l,^a;i/3ioga; / (l + 8M)-i/V2/3rf„ + 0(xV3(iogx)2/3+v^/6(loglogx)5/i2x 
Finally, by Lemma 14771 

o4/3 

Q(x) = ^-2-xi/Mogx + 0(a;^/=^(logx)2/3+^/6(loglogx)^/^'). 
31(3) 

□ 



5. Sums of multiplicative functions 

In this section we gather together the facts about sums of multiplicative functions that we used in 
Sections [3] and m All of the specific results we need are special cases of the following asymptotic 
formula, several variants of which have appeared in the literature. 

Proposition 5.1. Let g{n) be a nonnegative multiplicative function. Suppose that there exist real 
numbers U and n such that g{p'^) < U for all prime powers and 

> ^ = Klogw + Ogil) 

^-^ p 

p<w 

for all w > 2. Then the asymptotic formula 



n<y 

holds for all y > 2, where c{g) is the convergent product 

^Yf, , 9{p) , gip"^ 



r(K + i)-'-i\^ pj \ p p2 



Proof. This is exactly [6 , Proposition A. 3(a)] , except that we have replaced the hypothesis (n) -C 
for some constant « < 1/2" with "there exists U such that g{p°') < U for all prime powers p"'\ 
This is a strictly stronger hypothesis, however: any nonnegative multiplicative function satisfying 
g{p°') < U automatically satisfies g(n) < t/'^^") ^u,£ for every e > 0. □ 

In the following lemmas we use the standard notations 0(n) for the Euler phi-function, ii{n) for 
the Mobius mu-function, and uj{n) for the number of distinct prime divisors of n. Note that jj? is 
the characteristic function of the squarefree integers. 



Lemma 5.2. For all y >2,we have ^ \ / — — — ^ (log . 

n<y 



n0[n) 
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Proof. Define the multiplicative function 



9{n) 



(j){n) p-1'' 

p\n 



'n2'^(") -|--|- / 2p 



so that the sum we need to estimate is ^„<y g{n)/n. We note that g{p'^) = >/2p/ {p — 1) < 2 for 
all prime powers p", and we evaluate 



^ p ^ V » — 1 "p 

P<y ^ p<y * ^ ^ 

logp 



^yP^i-i/p 



p<j/ ^ p<y 

V2\ogy + 0{l). 



Proposition lS. 1 [ therefore applies with k = \/2, yielding 



E \h^^ = E ^ = ^^9)i}ogy)^' + 0,((logy)v^-i) « (logy)^ 
V nd)[n) ^-^ n 

n<y V ^ ' n<J/ 



as claimed (the ^-constant is absolute since the function g is fixed). □ 
Lemma 5.3. For ally > 2 and for any nonzero integer h, 



;3V JJV VP/ 



Proof. We sort the integers n <y according to their greatest common divisor d with h: 



n(j)(n) ^ ^ \ ndiin) ^ \ dm6(dm) 

n<y V ^ d\h n<y V ^ m<|//rf • ^ ^ 

{n,h)=d {m,h)=l 

The condition (m, /;,) = 1 implies that (m, c?) = 1, and so 



^ 2^(^){h,n) /2^W ^ / 2^M ^ /2^W ^ / 2'^(™) 



(m,?t)=l 
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The inner sum is ^ (log y)^ by Lemma [5^ and so it remains to bound the sum over d. But this 
is a multiplicative function of h, and so 



n (^+ v^^i^z^ 

^ * ^ i=0 ^ 



The lemma follows upon verifying that a/2/(]9 — 1)(1 — < 7/ ^ for all p > 2. □ 

For certain multiplicative functions g, we can actually find an asymptotic formula not just for 
J2m<y g{'m)/^ but also for X]m<?/ our next lemma we record only the upper bound, even 

though we can derive an asymptotic formula; in the lemma after that, the asymptotic formula is 
important enough to retain. 

Lemma 5.4. For all y >2,we have ^ H I ~^ ~F ) ^ ^" 

m<y p\m V / 

Proof. Define the multiplicative function g{m) = 7'^^'^^ jj,"^ (m) / y/m, where /i is the Mobius mu- 
function. One can check that 

p\m * d\m 

(because both sides are multiplicative functions of m, it suffices to check the equality on prime 
powers). We have 

E n (i + ^) = E Esw = E«M E 1 = !'E ^ + o(5;s(d)). 

m<y p\m y ' m<y d\m d<y d<y d<y 

d\m 

Since g{m) m^^/^+'^ for any £ > 0, the error term is Oe(?/^/^+^). We note that g{p°') < 7 for 
all prime powers p", and we evaluate 



g[.P) logP _ 71ogp ^ ^ 

p<y p<y 



Proposition 15.11 therefore applies with k = 0, yielding 



d<y 



The proposition now follows from equation (1271) and this last estimate. □ 
Lemma 5.5. Define S{y) = J2b<y Pi^)- For all y > 2, we have S{y) = -^ylogy + 0{y). 
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Proof. Define the multiplicative function g by its values on prime powers as follows: 
^b") = forpoddor a > 4; g{2) = 0, g{A) = 1, g{^) = 2. 

As in the proof of the previous lemma, one can check that p(m) = X]d|m ^^^^ 

s{y) = yY.'-^ + o(Y^3{d)). (28) 

d<y ^ d<y ^ 

Since (7 is bounded by 2, the error term is 0{y). We note that g{p°') < 2 for all prime powers p", 
and we evaluate 

V = V = logi/ + 1. 

p<y '2<p<y 
Proposition 15.11 therefore applies with k = 1, yielding 



= n (1 - ^) + Oil) = ^ + 0(1) = ^ log , + 0(1). 



The proposition now follows from equation (|28] ) and this last asymptotic formula. □ 

When p is generalized to p/ for any reducible, nonsquare quadratic polynomial / with integer 
coefficients, the proof of Lemma [53] generalizes to yield 

J2pfM=cif)y\ogy + Of{y) (29) 

m<y 

for some positive constant c(/). It is also easy to show (by splitting the sum into dyadic intervals, 
for example) that 

5:^«^"nogA (30) 

a<A 
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